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ABSTRACT 

This  thesis  is  concerned  with  the  steady  state  response  to  a  step 
input  of  a  non-linear  s ervomechani sm  having  viscous  friction  and  a  load 
position  unity  feedback  loop.  A  gear  train  with  "backlash  between  re- 
silient gear  teeth  was  located  between  the  motor  and  load. 

The  physical  equations  describing  the  system  are  presented  and 
adapted  to  a  form  for  the  Control  Data  Corporation  l60^-  digital  com- 
puter.  The  computer  program  used  in  this  analysis  is  presented  along 
with  its  flow  diagram  and  description  of  its  operation. 

The  results  of  this  thesis  are  presented  in  several  forms  using 
the  parameters  of  system  damping  coefficient,  distribution  of  friction 
and  inertia,  backlash  angle  and  coefficient  of  restitution.   The  re- 
sults are  applied  to  sample  problems  in  illustration  of  their  use  for 
design  and  analysis. 

The  authors  wish  to  express  their  appreciation  to  Dr.  George  J. 
Thaler  of  the  Department  of  Electrical  Engineering,  and  to  Dr.  William 
Wainwright  and  Mr.  Edward  N.  Ward  of  the  Department  of  Mathematics  for 
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1.   Introduction  and  Background 

The  problem  under  investigation  was  the  response  to  a  unit  step 
input  of  a  nonlinear  servomechanism  having  viscous  friction  and  load- 
position  unity  feedback  loop.  A  gear  train  with  backlash  between  the 
resilient  gear  teeth  was  located  between  the  motor  and  load.   A  de- 
termination of  the  possible  existence  of  a  limit  cycle  with  variation 
of  several  parameters  was  made.   The  variable  parameters  considered 
were  motor,  load,  and  system  time  constants,  backlash  and  the  co- 
efficient of  restitution  of  the  gear  teeth.   If  a  limit  cycle  existed, 
a  further  investigation  was  made  into  the  change  of  the  size  of  the 
limit  cycle  with  variation  of  the  above  parameters. 

If  a  gear  train  is  required  between  the  motor  and  load  in  a 
servomechanism,  the  gear  train  may  be  treated  as  ideal  with  perfect 
meshing  of  the  gears.   In  this  case  the  components  of  the  system  are 
joined  at  all  times  and  the  system  may  be  described  by  a  single  linear 
differential  equation.   The  classical  linear  solution  in  either  the 
frequency  or  time  domain  is  the  result. 

Present  production  methods  cannot  meet  the  requirements  of  an  ideal 
gear  train,  having  perfect  engagement  of  the  teeth,  nor  is  the  ideal 
gear  train  desirable  from  the  standpoint  of  wear  on  the  gear  faces. 
The  result  is  that  the  practical  gear  train  has  separation  of  the  gear 
teeth  when  the  velocities  of  the  motor  and  load  are  different.   The  re- 
sponse of  a  servomechanism  with  backlash  is  in  effect  discontinuous  or 
nonlinear.   The  net  system  acts  as  three  individual  but  dependent  systems 
with  the  boundary  conditions  of  each  being  specified  by  certain  physical 
laws . 

Two  accepted  methods  of  analysis  are  available  for  solution  of  the 


problem.  The  first  is  the  solution  by  describing  function  methods 
while  the  second  is  by  analytical  methods.   Describing  function 
methods  are  primarily  frequency  response  techniques  and  depend  upon 
the  assumption  that  the  nonlinear  element  may  be  considered  as  linear 
over  the  range  of  consideration.  A  further  assumption  is  made  that 
the  input  to  the  nonlinear  element  is  a  pure  sinusoid  with  no  harmonic 
frequencies.   The  describing  function  method  is  primarily  concerned 
with  steady  state  results. 

The  analytical  solution,  which  is  the  second  method  of  considering 
the  problem  is  primarily  a  ^transient  response  technique  concerned  with 
solution  in  the  time  domain.   This  approach  also  require.1-  certain 
assumptions  concerning  the  behavior  of  the  system.   The  assumptions 
used  in  several  investigations  are  pointed  out  in  the  following  para- 
graphs, while  the  assumptions  of  this  investigation  are  outlined  in 
Section  3>  Development  of  Equations.   The  validity  of  the  analysis  by 
any  method  depends  on  the  accuracy  with  which  the  real  system  has  been 
defined. 

Chestnut  and  Mayer,  Ref.  a,  have  considered  the  backlash  problem 
using  describing  functions  and  supported  the  findings  by  analog  com- 
puter studies.   The  main  difference  between  this  cited  work  and  previous 
describing  function  analyses  was  the  simulation  of  backlash  by  a  dead 
zone  between  the  motor  and  load  and  springiness  in  the  interconnecting 
shaft . 

Various  analytical  approaches  have  been  made  to  the  problem  of 
backlash,  Lutkenhouse,  Ref.  b,  used  graphical  methods  in  the  phase 
plane  to  solve  several  cases  of  plastic  impact  between  the  gear  teeth 
of  a  second  order  servomechanism.   Pastel  and  Thaler,  Ref.  c,  developed 


analytical  equations  for  the  existence  of  limit  cycles  using  phase 
plane  equations  for  the  case  of  plastic  impact  and  no  load  inertia. 
Knoll  and  Narud,  Ref .  d,  investigated  limit  cycles  in  the  phase  plane 
using  an  analog  computer.   The  investigation  of  Knoll  and  Narud 
covered  a  wide  range  of  parameters  for  the  case  of  plastic  contact 
between  the  gear  teeth. 

To  analytically  describe  the  total  system  response  of  a  second 
order  servo  with  "backlash,  three  differential  equations  and  one  or 
more  algebraic  equations  are  required.   The  system  as  a  whole  may  be 
treated  as  being  piece  wise  linear.   One  differential  equation  is  used 
to  describe  the  entire  system  when  the  gear  teeth  are  in  contact  and 
the  motor  is  driving  or  braking.   When  the  gear  teeth  are  not  in  con- 
tact, two  more  differential  equations  are  required,  one  for  the  motor 
alone  and  one  for  the  load  alone.   The  boundary  conditions  of  the  dif- 
ferential equations  are  determined  from  the  solution  of  one  or  more 
algebraic  equations,  expressing  the  laws  of  conservation  of  momentum 
and  energy.   For  plastic  impact  with  no  bounce  of  the  gear  teeth,  the 
law  of  conservation  of  momentum  must  be  satisfied.   For  perfect  elastic 
contact  between  the  gear  teeth  the  law  of  conservation  of  momentum  and 
the  law  of  conservation  of  energy  must  be  satisfied  simultaniously. 
Intermediate  cases  between  perfect  elastic  and  perfect  plastic  contact 
can  fulfill  only  the  law  of  conservation  of  momentum.   However,  total 
energy  accounting  may  be  made  in  the  intermediate  and  perfect  plastic 
impact  cases. 

New,  Ref.  e,  adapted  the  differential  equations  of  system  motion 
and  the  law  of  conservation  of  momentum  (plastic  contact)  to  a  digital 
computer  analysis  of  a  second  order  servomechanism  with  backlash.   The 


digital  computer  was  chosen  as  the  method  of  solution  for  this  investi- 
gation primarily  "because  of  the  flexibility  of  data  presentation  and  the 
ease  with  which  cyclic  parameter  variation  could  be  obtained. 

In  later  sections  of  this  thesis,  the  differential  equations  for 
the  system  and  the  algebraic  equations  for  the  impact  boundary  conditions 
are  developed.   The  computer  program  and  its  flow  diagram  and  modes  of 
operation  are  pointed  out.   The  results  of  the  solutions  of  problems 
are  presented  in  several  forms  and  observations  are  made  in  Section  5> 
Results  and  Discussion.   Sample  applications  of  the  results  to  mechani- 
cal systems  are  presented  fn  Section  6,  Application  of  Results. 

At  this  writing,  associated  work  in  the  area  of  transient  response, 
primarily  peak  c  ershoot  and  settling  time,  of  a  second  order  servo- 
mechanism  with  backlash  and  plastic  and  elastic  impact  is  being  pre- 
pared as  a  thesis  by  C.  E.  Andrews  and  R.  A.  Kelley  at  the  U.  S.  Naval 
Postgraduate  School. 


2.   Development  of  Equations 
The  equations  developed  were  those  of  a  phase  plane  analysis. 
Prior  to  defining  the  net  system  equations ,  the  assumptions  under  which 
the  analysis  was  made  will  "be  stated.   When  necessary  these  assumptions 
will  be  amplified  and  referred  to  later  in  this  work. 
It  was  assumed  that: 

1.  The  gear  teeth  were  initially  in  contact  and  the  initial 
conditions  of  the  system  were  all  equal  to  zero.   This  was 
later  proved  to  be  an  unnecessary  limitation  for  the  study 
of  steady  state  response. 

2.  Plastic  deformation  of  the  gear  teeth  during  steady  contact 
and  impact  and  any  torsional  deformations  of  driving  shafts 
are  negligible. 

3.  The  inertias  of  the  gears  and  drive  shafts  are  considered 
as  part  of  the  load  or  motor  inertia  depending  on  their 
attachment  in  the  system. 

4.  The  law  of  conservation  of  energy  was  completely  satisfied 
in  only  the  perfect  elastic  case  by  maintaining  the  total 
mechanical  rotational  energy  of  the  system  constant  at  the 
instant  prior  to  and  after  impact.   The  law  of  conservation 

of  momentum  is  satisfied  in  plastic ,    elastic,  and  intermediate 
cases.   When  the  law  of  conservation  of  momentum  is  the  only 
equation  required  to  be  satisfied,  the  energy  lost  from  the 
system  is  dissipated  in  the  heat  of  infinitesimal  deformations 
of  the  gear  teeth. 

5.  The  gear  teeth  are  in  contact  only  instantaneously  during 
impact  for  the  elastic  and  intermediate  cases  (excluding 


plastic)  and  that  the  Impulse  torques  of  drive,  friction  and 
"bearing  supports,  etc.,  are  zero  during  impact,  e.g. 


avAd      ©-=.   6- 


J  Tm  c\t  -  6 

The  coefficient  of  restitution  of  the  two  opposing  gear 


teeth  is  the  same  or  is  described  by  an  equivalent  co- 
efficient if  the  two  gear  teeth  are  of  unequal  coefficients. 

7«   Backlash  is  assumed  to  be  equal  at  all  points  on  the  gear 
circumference.   Backlash  is  measured  at  the  output  shaft. 

A  "block  diagram  of  the ''system  considered  is  presented  below. 
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The  equations  for  the  motor  and  load  in  contact  are: 
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The  equation  for  load  alone  with  gears  not  in  contact  is: 

Jl  6c+fL  6c  =  O 
or     ©c  +•  J^_  6C     =  O 

The  equation  for  the  motor  alone  with  the  gears  not  in  contact  is. 

The  method  of  analysis  used  by  New,  Ref .  e,  is  used  in  order  to 
be  able  to  examine  the  system  independently  of  its   MJ  n 
Defining  first      00  f\   T  =■  t 
and  differentiating      ^O^  UL  =  &T 

7F       dt*   dt-        dt* 
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at*-   (dt*ndtj    (j^)1 


•  * 


aZ^c      -     ©c*        Awcl       &,=  3/^M 

1  di-t 

C  4cCR    -    C  6c  di*  or  6c=  fi£ 

Using  the  equations  of  the  combined  system  and  the  load  alone, 

and  making  the  indicated  substitution  to  a  transformed  (*)  coordinate 
system  for  the  system  equation 

CO*  6*  +  a  t  Uii  &*  .+■  U)w  3<l  =  fon  £* 

ec  ^s  s*  +  '.at  -  £r 

and  introducting  the  slope  equations  of  the  phase  plane 

/&  4^  /  5§fi|?   t  fi&   =  6g   results 

If  the  load  equation  is  first  put  in  the  form  for  the  phase  plane 

and  then  transformed  to  the  (*)  coordinate  system  where  q  ^-  Q 

/yl     "&£,      ^  TL   "O  a        ~      O  results . 


The  same  result  may  be  obtained  by  setting  uj=\.      The  results  of  such 
a  transformation  require  inverse  scaling  for  practical  application, 
examples  of  which  given  in  Section  6.   It  is  noted  that  the  slope  of 
the  load-free  equation  is  mathematically  the  same  whether  in  the  phase 
plane  or  in  the  transformed  phase  plane.   It  is  pointed  out  at  this 
time  that  the  equations  later  developed  to  satisfay  the  laws  of  con- 
servation of  momentum  and  energy  are  independent  of  the  system  natural 
frequency. 

To  establish  equations  for  the  law  of  the  conservation  of  momentum 
and  a  relationship  satisfying  the  law  of  the  conservation  of  energy, 
the  following  schematic  is  used. 
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By  representing  the  inertia  of  the  load  and  motor  as  inertia  of  the 
gears  the  following  figure  is  obtained  at  the  instant  of  impact. 
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The  torque  equation  may  "be  expressed     '  =  0  u   —   J  y__0 
Impulse  is  equal  to  the  time  rate  of  change  of  momentum.   The  ex- 
pressions for  the  rates  of  change  of  momentum  of  the  gears  treated 
separately  may  be  written: 


**     — - 


m 


t 


av\A 


^VucH    -  ? w\ 


|  Fall  ^  (H    = 


6  m     i9 


w 


ToriKe  wart  or 


£ 


O, 


3,  ^d-t 

at 


ii 


where  the  primed  velocities  are  those  following  impact. 

If  the  momentum  is  to  be  conserved  in  the  system  the  impulse 
functions  are  equal  and  opposite . 

T™   j_  (  ^  -  d^)  =    f  *„  -    4c) 

from  substitution  of  the  equations 

The  assumption  that  all  other  impulse  functions  are  zero  at  the  instant 
of  impact  is  restated  at  this  point. 

A  definition  of  the  coefficient  of  restitution  will  be  made  for 
purposes  of  this  work: 

«.  *-    Urn  +  4'^        ^    _    (    p4l~  4! 


(4W  +  M  Cp4M-  ij 


It  will  be  shown  that  this  definition  of  e  will  satisfy  the  law 
of  conservation  of  energy.   Consider  the  case  of  e=0,     f 

in  which  case  for  plastic  impact  (e=0)  the  gears  are  moving  at  the 
same  velocity  following  impact. 

The  expression  for  the  coefficient  of  restitution  is  examined  for 
the  case  of  perfect  elastic  contact,  e=l,  for  which  the  law  of  conser- 
vation of  energy  is  satisfied. 

Solution  of  the  restitution  equation  for  e=l  yields 
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For  conservation  of  energy 
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The  momentum  equation  which  is  independent  of  e  is  \ 
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The  factored  energy  equation  is  now  divided  by  the  momentum 
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which  checks  with  the  definition  of  the  coefficient  of  restitution 

for  e=l. 

The  momentum  equation  implied  from  the  impulse  approach  that 

duT-^  0  hence        &c=    "O^   flftd  "d^5  6^.      The  definition  of  the 

coefficient  of  restitution  is  independent  of  position.   It  may  "be 

noted  that  the  equations  for  energy  and  momentum  depend  only  on  the 

inertia  distribution,  gear  ratio  and  angular  velocities. 

In  manipulation  of  these  equations  it  is  seen  that  load  inertia 

may  be  transferred  to  the7|  ,  shaft  by  multiplication  (~*     ,    the  same 

transfer  may  be  accomplished  with  load  friction.   Since  A  ,  and  0"  , 

are  related  to  U  \_      an(^   ^i   ^  may  ^e  reasoned  that  the  results 
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for  P  as  1  may  be  extrapolated  to  physical  systems  where  [^  *=  1. 

Since  a  total  energy  accounting  may  be  made  for  the  system  with 
the  equations  used,  the  rotational  energy  lost  is  attributed  to  the 
heat  of  deformation  of  the  gear  teeth.   This  transfer  of  energy  could 
in  fact  be  determined  from  the  equations  used.   Thus  the  conservation 
of  rotational  energy  e=l  is  a  special  case  of  a  broad  interpretation 
of  the  law  of  conservation  of  energy. 
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3.   Computer  Program  Development 
The  physical  equations  of  the  net  systems  were  programmed  for 
solution  using  the  Control  Data  Corporation  l66k   high  speed  digital 
computer  utilizing  paper  tape  program  input  and  magnetic  tape  output. 
An  IBM  717  line  printer  was  used  to  extract  data  from  the  magnetic 
tape.   The  Control  Data  Corporation  machine  library  and  the  U.  S. 
Naval  Postgraduate  School  computer  subroutine  library  were  used  for 
assembly,  Runge-Kutta-Gill  numerical  integration,  and  decimal  output. 
Several  changes  were  made  in  the  forms  of  the  physical  equations 
of  the  net  system  in  order  to  eliminate  duplicate  computing  operations 
and  to  fit  the  equations  to  a  form  suited  to  the  variable  parameters. 
In  the  table  of  symbols,  computer  mnemonic  (m  )  terms  which  are  used 
in  the  assembly  subroutine  have  been  indicated  by  parenthises  and  will 
be  defined  when  encountered. 

The  equation  of  the  system  with  motor  and  load  combined 

was  put  in  the  form 

6^       U)r\  £^>  -  6j*  ^c   -ft  6c      vuUere 

&  =      Z%  60n        a  Ad 

THETAL      =         Q^  ONE      =  Q)  fc 

4 

THETALD  m        \J(     -      Q^  ZETA   «*     § 

UDOT  =    Q^  OMEGANSQ  =    U) /\~ 

were  the  computer  m  terms  used.   Underlining  will  be  used  in  this 

section  to  denote  the  transition  to  the  m  terms  of  the  computer 

equations. 

The  differential  equations  were  solved  using  the  Runge-Kutta-Gill, 
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(RUNGS),  numerical  integration  method.   The  RUNGE  subroutine  re- 
quired the  definition  of  synonymous  terms  for  the  four  iterative 
cycles  used  to  produce  one  extrapolated  set  of  variables.   The  incre- 
ment of  the  independent  time  variable  chosen  was  0.01  sec.  for  all 
^}   QA     and  0.004   sec.  for  all  £  <^   Q,|  . 

The  load  free  equation 

f 
\3(^  +~  T*  L    0^-0  was  put  in  the  form 

C>l    i   o    ^    —  O  for  the  computer. 

A  separate  set  of  RUNGE  synonyms  was  used  for  the  load  free  equation 

in  order  not  to  destroy  the  previous  computations  before  they  were 

determined  of  no  further  value  in  computing  and  also  in  the  making 

of  program  decisions. 

Q  was  designated  V  and   ^-  'was   designated  VDOT  for  the  computer 

program. 

The  equation  for  the  motor  alone 

3"^  Tm  Tk/i 

was  put  in  the  form 

and  for  RUNGE,  the  terms 

*) 

Q  KVy  =  THETAM  -&u.  =  THETADM  =  W 

and      ^    =  WDOT       were  used. 

The  equation  representing  the  law  of  conservation  of  momentum  and 
the  definition  of  coefficient  of  restitution  were  combined  to  the  forms 
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for  the  computer  program.   No  provision  for  additional  m  terms  was 
made  to  denote  the  primed  values.   The  additional  terms: 

"were  used. 

Two  equations  were  used  to  define  the  "boundaries  of  operation  of 
the  load  and  motor  when  they  were  acting  separately, 

Two  major  decisions  of  the  computing  cycle  were  the  determination  of 
the  point  where  the  motor  and  load  would  float  free,  FLQATEST,  and 
the  response  of  the  system  to  impact  of  the  load  and  motor  gears  when 
the  boundary  conditions  on  position  were  met,  COMBTEST. 

The  first  decision,  the  point  of  float  free,  is  made  on  the 
equality  or  inequality  of  slopes  in  the  phase  plane,  i.e.,  at  separation 
TL=  H«.   Equating  the  slopes: 

The  above  equation  is  solved  after  computing  each  point  in  the 
combined  phase  trajectory.   If    /r<    '    fyj     /  0    "t^le  system  re- 
mains combined.   If   /\l <    /•  /\/>    "C*  Q    the  system  separates,  having 
the  response  motor  alone  and  load  alone  immediately  thereafter. 

The  second  major  decision,  that  of  whether  after  impact  the 
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response  should  "be  that  of  the  combined  system  or  the  response  of  the 
load  and  motor  acting  separately,  is  made  on  the  "basis  of  the  resultant 
velocities  after  impact  and  upon  which  side  of  the  backlash  the  motor 
and  load  positions  are  found. 

The  computer  program,  its  flow  diagram,  a  description  of  the  pro- 
gram operation,  and  an  explanation  of  the  readout  of  information  with 
examples  are  presented  in  the  Appendix. 
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k.     Methods  and  Scope  of  Examination 
Two  modes  of  operation  of  the  computer  program  were  used.   The 
first  mode  considered  was  the  phase  trajectory.   This  mode  was  used 
to  check  the  operation  of  the  computing  cycle.  Time,  velocity  and 
position  of  the  load  were  printed  when  the  system  was  combined.   When 
the  system  was  separated,  the  values  of  motor  velocity  and  motor 
position  were  also  printed.  Although  the  time  increment  used  for 
computations  was  either  0.01  sec.  or  0.00^  sec,  printed  outputs 
for  the  first  mode  were  taken  only  every  0.1  sec.  Additional  printed 
outputs  were  taken  at  the  time  of  contact  of  the  gears  just  prior  to 
and  immediately  following  the  solution  of  the  momentum  and  restitution 
equations.  A  typical  printout  of  the  phase  trajectory  mode  of  operation 
is  shown  in  Appendix  C. 

The  second  mode  of  operation  was  used  for  printing  only  the  maxi- 
mum computed  positive  load  overshoot  position  for  each  cycle,  the 
associated  load  velocity  and  the  exact  problem  time  of  the  computation. 
A  sufficient  number  of  print-outs  for  each  problem  was  obtained  to 
prove  limit  cycle  existence  and  average  size.   Since  this  investigation 
was  concerned  primarily  with  the  steady  state  response,  the  second  mode 
of  operation  was  the  one  utilized  to  obtain  the  majority  of  the  data. 
This  read-out  method  markedly  decreased  the  data  reduction  time  for  the 
problem  analysis,  as  opposed  to  the  analysis  of  steady  state  conditions 
provided  by  a  full  phase  trajectory  print-out.   Since  the  time  required 
for  computer  read-out  was  several  magnitudes  greater  than  the  computing 
time,  valuable  computer  time  was  saved  by  this  mode  of  operation.   An 
example  of  the  printed  output  obtained  when  only  maximum  load  position 
overshoot  was  of  interest  is  given  in  Appendix  0.     Approximately  1^00 


19 


phase  trajectories  were  solved  using  this  mode  of  operation.   Each 
solution  required  an  average  of  three  minutes  of  computer  time. 

For  the  purpose  of  general  examination  of  the  problem  it  was 
assumed  that  P  =  1,  f^  =  1,  6J^  =  1,  3^  -P  3L  =  1.   Specific  solutions 
were  made  with  0  £  1  }   (j(L^  1  to  determine  scaling  effects.   The  re- 
sults of  these  solutions  are  presented  in  Section  6,  Application  of 
Results . 

The  program  can  also  he  made  to  solve  linear  systems  by  setting 
A=   0. 

The  major  variable  parameters  used  in  this  investigation  were: 

<J  m   0.1,  0.2,  0.3,  O.k,    0.5,  0.6,  0.8,  1.0 
(°\=  0   ,  0.2,  O.k,    0.6,  0.8,  1.0 

VJ^=  o.i,  0.2,  0.5,  o.8,  0.9 
■3^  0.9,  "oTB7,  0.5,  0.2,  0.1 

Z\  =0.3  radians 

e  =  0,  0.6,  0.8,  1.0 
Additional  parameters  of 

\  \ l  =  0.0U,  0.1,  0.9;  and   A  «=  0.01,  0.02,  0.01,  0.15  radians 

Ft 

were  used  in  certain  instances  to  examine  particular  characteristics  of 
the  response. 

The  program  was  limited  somewhat  in  that  the  value  of   P _£jr  =  ^q 
and  zero  were  exluded  due  to  generation  of  undefined  mathematical 
quantities. 

The  value  of  backlash  was  made  abnormally  large,  A  =0.3  radians, 
to  allow  an  easier  interpretation  of  the  non-linear  response.   Since 
the  influence  of  J^    was  linear,  this  caused  no  inaccuracies.  All 
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graphs  for  limit  cycle  size  were  plotted  with  /\     =  0.3  radians. 
The  effect  of  values  for  ^\    other  than  0.3  is  also  shown  in  graphi- 
cal form. 


21 


5.  Results  and  Discussion 

Two  types  of  system  response  were  observed.   The  first  type  was 
a  convergence  to  zero  error,  which  was  a  case  of  static  and  dynamic 
equilibrium  of  the  system.   This  is  a  characteristic  of  a  linear 
second  order  servomechanism.   The  second  type  of  system  response 
observed  was  the  divergence  or  convergence  to  a  state  of  dynamic  equi- 
librium.  This  second  type  of  response  is  termed  a  limit  cycle  and  is 
characterized  by  cyclic  travel  of  the  load  through  the  same  points 
of  the  phase  trajectory.   When  the  second  mode  of  computer  operation 
was  used,  both  the  maximum  load  position  per  cycle  and  the  period 
of  the  recorded  position  were  analyzed  to  determine  if  the  character- 
istic response  was  a  limit  cycle.   Exact  repetition  of  cyclic  values 
was  impossible  due  to  the  numerical  methods  used  and  the  computer 
round-off  error  for  the  output  routine.   Typical  examples  of  the 
limit  cycle  and  no  limit  cycle  response  are  given  in  Appendices 
D  1  and  D  2. 

No  case  of  dynamic  instability  or  divergence  of  the  system  with- 
out bound  was  observed  in  the  investigation.   Intuition  might  lead 
to  this  same  conclusion  when  the  following  unique  characteristics  of 
this  system  are  pointed  out:   (a)  As  a  maximum  limit,  energy  was 
conserved  at  the  instant  of  impact.   (b)  Energy  was  supplied  to  the 
system  only  in  finite  amounts  and  at  a  finite  rate  by  the  error 
detector.   (c)  Energy  sinks  were  present  in  the  system  frictional 
elements  while  energy  storage  units  were  present  in  both  the  motor 
and  load.   (d)  No  undefined  limits  appeared  in  the  applied  physical 
equations, 


P  jjL-  =   (2>  and  qQ  excluded. 
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The  results  of  the  investigation  are  offered  in  the  form  of 
charts  for  parameter  areas  of  limit  cycle  existence  and  nonexistence 

in  Fig.  1  for  e  =  0. ,  0.6,  0.8  and  Fig.  2  for  e  =  1.0.   The  para- 

"2  r  2. 

meters  P    \ ,_  and   D  Tl  are  designated  abscissa  and  ordinate 

respectively.   The  parameter  points  of  examination  for  which  solutions 
were  obtained  are  circled.   The  zone  enclosed  by  a  particular  jf  line 
is  largely  an  interpretation  by  the  authors  of  the  results  of  the  in- 
vestigation, and  shows  the  area  in  which  no  limit  cycle  existed.   Dae 
to  the  large  change  of  the  variable  parameter  values  between  solutions, 
the  charts  are  not  exact.  All  circled  intersections  interior  to  the 
\    line  indicate  parameters  which  resulted  in  no  limit  cycles,  e.g., 

Fig.  1  for  e  =  0,  0.6,  0.8  at  a   ^   =  O.k,    the  values  of   ^  k   = 

FT 

0.6,  0.8  at  a  P  ^3j_  =  9.0  and  the  value  of   ^  Vl.   =  0.6  at  a 

~Zyv\  ^T 

0  PL  =4.0  resulted  in  no  limit  cycle.   Similarly,  for   "\  =0.5, 

the  values  of  £  V  i_  =  O.k,   0.6,  and  a   O  Jl  =  1.0  as  well  as  all 

other  parameter  intersections  interior  to  the   >T  =0.5  curves,  re- 
sulted in  no  limit  cycles. 

All  circled  parameter  intersections  exterior  to  particular   \  lines 
define  paramters  which  were  examined  and  resulted  in  a  steady  state 
limit  cycle.   The  area  to  the  right  of  the  single  \     =1.0  line 
designates  the  parameters  for  which  no  limit  cycles  were  observed. 
Limit  cycles  were  observed  for  all  A  <C    0.3  for  the  parameters  examined. 

The  investigation  was  made  with   ^  =0.3  radians;  however,  re- 
sults of  the  investigation  of  variation  of  limit  cycle  size  indicate 
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that  existence  of  the  limit  cycle  is  independent  of  the  size  of  back- 
lash in  a  nonlinear  system  for  a  A  greater  than  zero.   Variation  of 
/\   had  the  effect  of  a  linear  variation  of  limit  cycle  size  only. 

The  zones  depicted  on  the  existence  charts  are  approximately 
represented  by  straight  lines;  however  there  is  reason  to  believe  that 
a  complete  examination  would  show  a  curvature  in  the  lines.   This 
observation  will  be  elaborated  upon  later  when  the  relationship  between 
the  existence  charts  and  the  figures  of  limit  cycle  size  is  discussed. 

From  Fig.  1  it  is  seen  that  there  is  only  a  very  small  difference 
in  the  limit  cycle  existence  zones  for  plastic  impact  (e  =  0)  and  the 
intermediate  values  of  the  coefficient  of  restitution,  e  =  0.6  and  0.8. 
When  comparing  the  elastic  impact  (e  =  l)  in  Fig.  2  to  the  plastic  and 
intermediate  cases  in  Fig.  1,  the  most  obvious  point  of  comparison  is 
the  drastic  reduction  of  the  areas  in  which  limit  cycles  did  not  result 
for  e  =  1.0.   However,  it  may  be  generally  observed  that  with  any  e 
held  constant,  the  area  in  which  a  limit  cycle  could  not  be  obtained 
increased  with  the  increase  of  system  \T     . 

For  all  values  of  e,  there  is  an  apparent  common  enclosed  zone 

centroid  which  is  described  by  the  system  parameters  ( J  f  /F  =0.5 

and  P  J  /J  =1.0.   These  parameters  result  in  the  system  time  con- 
\   L/  m 

stants  fT  /J.  =  f  /  Jm,  since  Cx)      =  1  and  P  =  1  were  chosen  for  the 
1/  L    a/    '        wn         \ 

investigation.   There  is  also  an  apparent  symmetry  about  the  centroid 
of  opposite  quadrants  which  can  be  related  to  the  fact  that  since 
60  =  1  and   D  =  1,  the  points  diametrically  opposite  across  the 
centroid  of  the  figure  (using  the  decimal  values  of  the  parameter  inter- 
section) represent  an  exact  exchange  of  component  time  constants.   (See 
Fig.  2  for  example).   Two  diametrically  opposite  points  may  have 
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exactly  the  same  motor  and  load  time  constants;  however,  they  are 
descriptive  of  two  different  systems  since  the  friction  and  inertia 
distributions  are  different  for  the  two  points. 

Figs.  3  through  26  present  the  maximum  positive  error  in  radians 
of  the  limit  cycle  for  a  unit  step  input  and  with  backlash  equal  to 
0.3  radians.   The  use  of  the  charts  is  explained  in  Section  6,  Appli- 

cation  of  Results.   Each  figure  is  presented  with  p  f T /F  and  e 

1 

constant.   The  independent  parameter  Q  J  ,J  is  indicated  as  abscissa 

\  L/m 

and  curves  of  constant  3  are  plotted  with  limit  cycle  error  as  the 
ordinate.   The  results  were  presented  on  a  log  log  plot  since  the 
observed  limit  cycle  errors  encompassed  a  range  from  0.01  to  1.0 
radians,  and  the  appearance  and  disappearance  of  limit  cycles  and 
the  change  of  limit  cycle  error  was  more  readily  apparent  on  the  log 
log  form. 

In  general  as  e  increases,  the  size  of  limit  cycle  decreases,  all 
other  parameters  remaining  constant.   However  this  is  not  without  a 
few  exceptions -the  generalization  is  most  accurate  for  large  size  limit 

cycles.   It  may  be  noted  that  when  viewing  the  individual  charts,  if  a 

7. 

limit  cycle  exists  at  very  low   P  J  /j  ratio  there  is  a  tendency  for 

\   L  M 

the  same  system  to  have  a  smaller  limit  cycle  for  higher  P  j/j  ratio. 

\  L  m 

In  addition,  for  these  same  conditions,  if  limit  cycles  exist  for 
several  values  of  ^  ,  there  is  a  convergence  toward  the  same  size  of 
limit  cycles  for  each  \  as  Q>  J  /j  increases.  At  the  value 

*  /Fm 
i/  T 

5- 


D  fT/Fm  =  0,  a  limit  cycle  will  always  occur  for  any  value  of  e  or 


"2. 
Viewing  the  charts  in  sequence,  with  either  e  constant  or  Of  /F 

constant,  it  is  noted  that  there  is  continuity  of  pattern  flow  between 
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adjacent  charts.   This  continuity  "between  charts  appears  to  form  a 

symmetrical  pattern  of  constant  "x  lines  which  can  he  interpolated  for 

2. 
intermediate  values.   When  the  figures  for   O   f  /F  =  0.^   and  0.6 

\   L  t 

are  analyzed,  this  symmetrical  pattern  is  observed  in  each  chart  about 

*-   . 

PJ  / J  =  1.   The  limit  cycle  decreases  and  disappears  at  the 
L  m  _ 

opposite  extremes  of   p  J  /J  .   The  charts  of   Df  /F  =  0,k   and  0.6 

1   L  m  \  L  T 

"2_ 
are  symmetrical   to  each  other,  while  P  fT/FT  =0.2  and  0.8  are 

similar.   At  the  extreme  values  of  fa  f  /F  =  0  and  1.0,  no  symmetrical 


comparison  can  be  made.  An  additional  point  of  symmetry  about 
ft  fT/F  ■  0.5  can  be  visualized  for  all  values  of  e.   These  para- 
meters for  symmetry  are  noted  in  the  discussion  of  the  existence  charts, 


The  irregularities  in  the  curves  of  medium  values  of  *N  can  not  be 
explained  by  the  authors;  however,  it  may  be  mentioned  again  that  the 
points  of  symmetry  are  characterized  by  exact  exchange  of  component 
time  constants. 

The  noted  symmetry  of  the  limit  cycle  size  charts  shows  a  direct 
relationship  to  the  symmetry  found  in  the  existence  charts.   Since  the 
existence  charts  show  contours  of  constant  ^r  for  a  surface  of  zero 
limit  cycle  size,  one  is  led  to  the  possibility  of  a  three  dimensional 
presentation  to  show  a  variation  of  limit  cycle  size  for  the  same  para- 
meters used  in  the  existence  charts.   Such  a  three  dimensional  figure 
for  constant  e  using  D  f  /F  ,  £>  J  /j  and  j  as  the  axes  is  shown 
in  Fig.  27.   Other  surfaces  could  be  added  above  this  zero  limit  cycle 
"floor"  to  show  the  variation  of  a  given  limit  cycle  size  by  contours 
of  ^?  .   The  ultimate  result  of  this  line  of  investigation  would  be 
the  evaluation  of  the  poles  and  zeros  for  a  conformal  mapping  presen- 
tation.  Unfortunately,  time  did  not  permit  the  authors  to  investigate 
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this  avenue  further. 

The  results  of  the  investigation  of  the  dependence  of  limit  cycle 
size  on  the  amount  of  "backlash  are  presented  in  Figs.  28  and  29  using 
coefficients  of  restitution  of  0  and  1.0  at  a  \  =  0.6.  A  thorough 
investigation  was  made  at  this  \*    since  it  is  an  average  system 
parameter  often  used  in  design  work.   The  limit  cycle  size  was  found 
to  vary  directly  with  "backlash  size  without  exception  between  the 
values  of  0.3  and  0.01  radians.   When  limit  cycles  existed  at  £^  =  0.3, 
there  was  no  disappearance  of  limit  cycle  with  decreasing  magnitude  of 

.A  .   Where  no  limit  cycle  existed,  at  f\    =  0.3,  limit  cycles  con- 
tinued to  be  non-existent  throughout  the  range  of  _^\  examined. 
Several  solutions  of  problems  with  ^  =  0.6  and  e  =  0.6  &0.8  and  "V*  ■ 
0.1  &  e  =  1.0  &  0  with  various  /\^    values,  proved  the  linear  variation 
of  limit  cycle  size  with  backlash  size  for  all  values  of  e. 

When  a  serve-mechanism  under  steady  state  conditions  is  disturbed 
by  a  small  perturbation,  it  will  return  to  the  original  steady  state 
conditions.   However,  the  transient  response  due  to  the  perturbation 
may  be  quite  different  from  the  transient  response  which  occured  due  to 
the  original  signal  input.   The  concepts  of  energy  and  static  and  dy- 
namic equilibrium  mentioned  in  Section  1  may  be  better  understood  as  a 
result  of  the  perturbation  technique  if  the  usual  initial  conditions  of 

■Ov  (o)  =  cl  (o)  =  0  are  considered  as  a  disturbance  from  the  stability 
conditions  predetermined  by  selection  of  system  parameters  and  signal 
input. 

This  same  small  perturbation  technique  was  simulated  mathematically 
using  the  digital  computer  to  obtain  several  solutions.   The  cases 
examined  and  their  initial  conditions  were: 
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a)  ^  =  0.6,       Q  JL/Jm  =  0.8/0,2,        ft  fL/FT  =  0,    e  =  0,      A  =  0.3 

-6cc    (0)   =  0     -6  c    (0)   =  1.2 
4    m   (0)   =  0       &m   (0)   =  1.0 
system  separated. 

b)  ^J  =  0.6,     f  JL/Jm  =  0.1/0.9,         PfL/FT  =  i-0^  e  =  !•<>>    A    =  0.3 
^  c   (0)   =  1.0  £  c    (0)   =  1.0 

-6  m   (0)  =  1.0    ,    &   m  (0)  -s  1.0 
system  combined 

c)  *§  =  OA,      9  JL/Jm  ■  0.5/0.5,       PfL/FT  =  °'h>  e  =  °'6'    ^  =  °'5 
■^  c    (0)   =  0  ;     S  c   (0)   =  1.005 

$   m   (0)   =  0        £>m  =  1.00 

system  separated 

In  the  three  cases  examined,  the  same  limit  cycle  was  obtained, 
both  in  size  and  period,  as  that  resulting  from  the  normally  used 
initial  conditions,  ^-c  (o)  =  (j  c  (0)  =0.   The  maximum  difference 
between  the  perturbation  response  and  the  original  limit  cycle  size 
for  the  three  cases  was  18  x  10~y  radians.   Two  types  of  transient 
response  were  observed.   These  were  convergence  to  limit  cycle  in  the 
first  two  cases  and  divergence  to  limit  cycle  in  the  last  case. 

The  results  obtained  by  the  small  perturbation  technique  supported 
the  validity  of  the  analysis  of  the  problem  by  digital  methods.   This 
result  also  concurs  with  the  findings  for  the  case  of  e  =  0  by  Knoll 
and  Narud,  Ref .  d,  in  that  the  steady  state  limit  cycle  size  is  completely 
specified  by  the  choice  of  system  parameters  and  independent  of  system 
initial  conditions. 

Mention  is  made  here  that  the  computer  program  could  not  be  run 
normally  for  the  perturbation  method.   The  computer  was  stopped  after 
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the  print  out  of  system  constants,  then  the  initial  conditions  of  the 
servo  system  were  inserted  manually.   The  program  was  then  restarted 

at  the  desired  point  for  computing,  either  as  a  combined  system  or 

with  the  motor  and  load  separated.   The  program  was  then  allowed  to 

operate  to  the  normal  completion  of  the  solution. 

The  work  of  Ref .  d  reported  that  the  system  transient  response 
for  the  case  e  =  0  arrived  at  a  steady  state  limit  cycle  by  two 
possible  phase  trajectories.   One  of  these  phase  trajectories  was 
smooth  convergence  to  the  limit  cycle,  while  the  second  was  a  converging 
overshoot  toward  the  inside  of  the  limit  cycle  followed  by  a  divergence 
to  the  limit  cycle.   From  the  results  of  this  investigation  for  the 
cases  e  =  0,  0.6,  0.8  and  1.0,  similar  transient  indications  were  ob- 
tained. 

Several  computer  solutions  were  accomplished  for  the  cases  J^    =  0.3 
radians,  e  =  0,  1.0,   D  -  1.0,  /j)    n  =  1.0,  K  -  J  '  £   1.0.   In  addition 
various  combinations  of  O  f   1.0  and  ^  n  /  1.0  were  obtained.   The 
purpose  of  these  solutions  were: 

1.  To  determine  the  general  applicability  of  the  computer  program. 

2.  To  verify  the  scaling  developed  in  Section  2,  Development  of  Equations. 

3.  To  determine  the  applicability  of  the  results  of  this  thesis. 
The  following  results  were  determined: 

A.  Two  solutions  were  obtained  (e  =  0&  1.0)  for  the  case  ^  =  0.4, 
P-1.0  U)n   =  1.0,  J  =  2.0,  K  =  2.0,   (^£^=0.4,   PJl=4.0. 
The  results  of  these  solutions  were  the  same  (time,  load  velocity,  load 
position,  and  limit  cycle  size)  as  the  solutions  obtained  for  \JU    -   0.4, 

?3i-  =  4.0,  ^  =  0.4  using  the  nondimensional  parameters  [/)   n  =  1.0, 
P  -  1.0,  J  -  1.0,  K  =  1.0. 
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Bo   Solutions  were  obtained  for  e  =  0  &  1.0  for  the  case  "X  =  0.4, 

P  =  1.0,  00  n  -  2.0,  J  =  1.0,  K  =  4.0,   P  Lb  =  0.4,   5?Ti_  =  4.0. 

The  results  of  these  solutions  yielded  the  same  size  maximum  overshoots 

and  limit  cycle  size  as  the  nondimensional  solutions  for  ^  =  0.4, 

\  t L  =  0.4,    °  ."J  L  =  4.0.  The  times  of  maximum  overshoots  in 
\  Ft  ^Jy^ 

this  case  were  one=half  the  times  of  occurence  of  the  same  overshoots 

for  the  nondiniensional  solution.   The  velocity  at  a  maximum  overshoot 

was  twice  the  value  of  velocity  obtained  for  the  same  overshoot  in  the 

nondimensional  solution.   This  same  comparison  can  he  made  for  Cases  A 

and  B.   The  results  of  these  *and  other  solutions  prove  the  validity 

formulae  used  in  Section  2,  Development  of  Equations,  concerning  the 

transformation  to  the  coordinate  system.  (*) 

C.  Solutions  were  obtained  for  the  case  e  =  0  &  1.0,  ^3  =  0.4,  P  -  0.5> 

(j)n  =  1.0,     _£l  =  0.8,  Tl  =  4.0,  \[±   =  0.2,   P_3j_  =1.0.   The 

Ft       ~5*n  Ft        ~3">^ 

results  of  these  solutions  were  the  same  (time,  load  velocity,  maximum 

overshoots  and  limit  cycle  size)  as  for  the  nondimensional  solutions 

-^  =  0.4,  P  =  1.0,  &)n  =  1.0,   ^£l  =  0.2,    P_I^=  1-0. 
^  Ft  T>v\ 

D.  Solutions  were  obtained  for  the  case  e  =  0  &  1.0,  "^  =  0.4, 

P  -  0.5,  (A   n  =  2.0,   £t-  m   0.8,   3j-  =  4.0,    3-ji-    =0.2,  ^^Ju  =  1.0. 

PT         "T/va  p  T  "J  w\ 

The  results  of  these  solutions  yielded  the  same  size  maximum  overshoots 

and  limit  cycle  size  as  the  nondimensional  solutions  for  ~\     =0.4, 

pyL_=  0.2,  P  jvK  =  1.0.   The  times  and  velocities  in  this  case  were 


Ft   ,  ~Zw 

scaled  by  LA)  n  as  in  Case  B. 

Several  phase  trajectories  using  e  =  0  were  compared  with  those 
obtained  by  New,  Ref .  e.  Generally  good  correlation  was  obtained  on 
the  value  of  limit  cycle  error  and  in  the  transient  response.  Small 
disagreements  were  expected,  resulting  primarily  from  differences  in  the 
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programming  procedures  and  analysis  of  results  "by  read-out  methods. 

Sample  phase  trajectories  for  the  cases  of  e  =  0  &  1.0  are  pre- 
sented in  Figs.  30  and  31«  Although  the  transient  analysis  of  the 
system  is  beyond  the  scope  of  this  work,  these  figures  may  aid  in 
understanding  the  physical  problem. 

Of  academic  interest  was  the  accidental  operation  of  the  computer 
program  for  several  solutions  using  \    =  0.   Several  responses  were 
observed.   One  was  a  nearly  pure  oscillatory  system  from  the  first 
overshoot  and  thereafter  another  was  a  very  slow  convergence  or  slow 
divergence  of  the  oscillations  after  the  first  overshoot.   The  con- 
vergence or  divergence  probably  resulted  from  errors  generated  in 
numerical  iterations  and  approximations. 
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6.  Application  of  Results 

In  this  section,  the  results  of  this  thesis  will  be  applied  to 
hypothetical  physical  systems.  Methods  will  be  indicated  by  which  physical 
systems  may  be  analyzed  or  designed.  Attention  is  called  to  the  assump- 
tions  stated  in  Section  2  under  which  this  study  was  made. 

It  is  desirable  to  find  an  equation  relating  certain  system 
variables,  therefore  the  equations: 


(i) 


D 


A      - 


]S—  -       P  k 


a  v\  d 


(2) 


3.t0dA    -      ^  t|3  Pl 


? 


I^L 


p 


T  w     4-   I 


Tm+  PZ^i 


will  be  manipulated  to  the  form 


~]VA 


[ 


+  to 


A  w\ 


H 


(3) 


^  =   ±  feu    [P^u         J       VTv^       ^|  +  fl 


P   0L 


or 


fc) 


S- 1  Hfe 


r*        4.  | 


£ 


'/» 


In  view  of  the  form  of  Equation  (!»■)  for  \  and  recalling  that  the 
variable  parameters  used  on  the  figures  were   P  Ai:  and    \^  3J=.  = 


1 


in  fact  includes  the  other  two  variable  parameters 


it  is  seen  that  the  expression  for  *P 
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From  discussion  concerning  the  effects  system  natural  frequency, 
it  will  "be  recalled  that  the  value  of  overshoot  and  hence  limit  cycle 
are  not  affected  by  (jj    n>   however,  the  time  and  velocity  are  scaled 
for  (J{j   n  f   1.0.   Thus  the  transient  response  characteristics  may 
he  related  to  the  same  parameters  used  in  these  applications. 

If  a  motor  and  load  are  selected  for  a  system  and  the  limitation 
is  imposed  that  system  operation  shall  not  result  in  a  limit  cycle, 
Equation  (k)   and  the  limit  cycle  existence  charts  may  be  used  to 
determine  the  unspecified  system  parameters   "X  ,  K  and  \~     .   The 
coefficient  of  restitution  for  the  proposed  gear  train  material  may 
be  determined  by  the  use  of  the  equation  in  Section  2. 

To  proceed,  select  a  value  of  K  and  P  ,  then  determine  the  value 
of  sT^  from  Equation  (k).      Enter  the  limit  cycle  existence  chart  for 
the  value  of  e,  system   \  and  the  parameters  r  ^ L  and  -- — L . 

Determine  if  the  selected  parameters  describe  a  point  of  no  limit 
cycle.   Since  limit  cycle  existence  is  independent  of  backlash  (^  ^  0 
excluded),  backlash  is  not  a  variable  parameter-   Successive  trial 
values  of  the  unspecified  parameters  may  be  required  to  meet  the 
limitation  of  no  limit  cycle.   If  it  is  also  required  that  certain 
transient  response  characteristics  be  met,  these  conditions  could  be 
examined  by  a  similar  method. 

In  the  event  that  no  acceptable  parameters  are  found  to  satisfy 
the  system  requirements  for  no  limit  cycle,  the  figures  for  limit  cycle 
size  may  be  examined  to  obtain  a  minimum  size  limit  cycle.   The  method 
of  solution  is  similar  to  that  previously  described. 

Example:   The  material  intended  for  use  in  the  gears  has  an  e  =  0.6, 
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The  motor  and  load  have  the  parameters  fm  =  0=64,  Jm  =  0,25,  f  =  0.64, 
JL  =  1.0  and  it  desired  that  the  system  operate  with  no  limit  cycle. 
If  the  values  are  selected  K  =  4,0,  P  =  0.5,  the  values  p  "Ju  =1.0 
and  "W      =  4,0,   1 =  0.2  result. 

Solve  Equation  (k)   for  "^  ■  .5(.64)(5+l)      •  5 (  thus 

*~^   =  0.4,   Enter  the  chart  for  limit  cycle  existence,  Fig.  1, 


e  =  0.6,  ^  -  0.4,   —I =  0.2  }         fi  Jl. 


^  *'  3 


1.0 


It  can  "be  seen  that  limi-t  cycle  will  exist  for  these  parameters. 

Since  K  is  most  easily  varied  in  equation  (4),  a  new  K  is  selected 
at  1.0  in  an  effort  to  find  parameters  for  no  limit  cycle.   It  is  seen 
that 


K-  u  K'  H.o  V     i 

with  all  other  parameters  constant,  thus   ^  =0.8       Enter  Fig.  1 


with  e  =  0.6 

J 

These  parameters  describe  system  operation  with  no  limit  cycle  and  are 
well  within  the  area  limits.   The  transient  response  of  this  system 
may  be  undesirable  and  various  combinations  of  p  and  K  would  have  to 
be  tested  to  satisfy  the  additional  requirements  of  the  problem.   It 
will,  be  restated  that  the  existence  areas  of  Figs.  1  and  2  are  not  exact. 

The  figures  of  limit  cycle  size  may  be  used  in  the  same  manner  as 
the  existence  charts.   The  same  parameters  of  the  first  case  are 
assumed:   K  =  4.0,   P  =  0.5,    P  Ju  =  1.0,  _! =0.2   ^  =  0,4 
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e  =  Oo6  ^\   =  0.1  radians.   Enter  Fig.  10  with  the  parameters  given. 
The  ordinate,  Limit  Cycle  Error  /  0.3  \  =  0.13  radians,  is  obtained. 

V  A  / 

The  maximum  positive  error  in  the  limit  cycle  =  0.13  (O.l)  =  0.0^3^ 

>p  To3T  dV 

radians.   Enter  the  same  figure  for  e  =  0.6,  \  =0.8      (  Ji_  =  1.0  . 
■f^   l.  \       =   0,2  .  /\  =  0=5  radians;  this  is  the  second  case  of  exami- 
nation  for  existence.   No  \^  =  0.8  curve  exists;  hence  no  limit  cycle 
exists  for  this  choice  of  parameters  for  any  nominal  value  of  "backlash. 
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7.   Conclusions 
Within  the  limitations  of  the  assumptions  made  and  the  scope  of 
this  thesis,  it  is  concluded  that: 

A.  Two  types  of  system  steady  state  response  result.   The  type  of  re- 
sponse is  Independent  of  backlash  size,  system  natural  frequency  and 
initial  conditions.   The  types  of  response  are  convergence  to  signaled 
input  position  and  convergence  to  a  limit  cycle. 

B.  If  a  limit  cycle  exists  in  the  steady  state  response  of  a  system, 
the  size  of  the  limit  cycle  varies  directly  with  the  size  of  backlash, 
all  other  parameters  constants 

C.  The  greatest  change  in  system  steady  state  response  occurred  for 
the  case  when  system  rotational  energy  is  conserved  (e  =  1.0). 

D.  A  point  of  symmetry  for  the  type  of  system  steady  state  response 
exists  with  equal  motor  and  load  time  constants.  For  medium  to  high 
values  of  system  damping  coefficient,  symmetry  with  respect  to  limit 
cycle  size  exists  for  an  exchange  of  motor  and  load  time  constants. 

E.  The  results  of  this  thesis  are  presented  in  a  form  which  may  aid 
in  system  design  or  analysis. 
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8.   Becommendations 
From  the  results  of  this  thesis  it  can  be  seen  that  there  remain 
many  fertile  areas  for  further  investigation.   For  studies  of  this 
type  it  is  strongly  recommended  that  the  digital  computer  be  used 
because  of  its  speed  and  versatility.   The  limits  of  investigation 
depend  to  a  great  extent  on  the  ingenuity  of  the  programmer.   Some 
recommended  areas  of  future  investigation  are: 

a)  More  precise  definition  of  the  limit  cycle  existence  zones 

and  a  more  thorough  coverage  of  the  limit  cycle  size  of  similar  systems 
having  the  same  type  of  parameters.  A  three-dimensional  presentation 
as  mentioned  in  the  discussion  could  be  obtained. 

b )  Include  torsion  and  other  deformations  of  the  mechanical  system 
in  the  physical  equations  of  the  system  to  determine  their  affect  on 
the  system  response. 

c)  An  analysis  of  the  limit  cycle  existence  zone  and  transient 
response  of  a  system  to  determine  if  there  is  a  relationship  between 
the  two.  An  analytical  expression  for  the  possibility  of  limit  cycle 
existence  might  be  developed  in  terms  of  system  parameters, 

d)  Examine  the  transient  and  frequency  response  of  the  system  due 

to  various  Q(?  input  functions.   This  would  require  minor  modifications 
to  the  computation  section  of  the  existing  computer  program. 

e)  Since  the  physically  realizable  quantities  of  position,  velocity 
and  acceleration  of  both  the  motor  and  load  are  available  as  computed 
quantities,  the  study  of  the  compensation  problem  and  its  optimization 
are  obvious  areas  for  further  consideration. 

f )  The  general  study  of  other  higher  order  nonlinear  systems  is 
suitable  for  digital  computer  programming  since  most  numerical  inte- 
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gration  methods  will  handle  systems  having  equations  of  any  order. 
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Maximum  Error  of  Limit  Cycle  from  Unit  Step  Input 
Limit  Cycle  Error  (Radians)  \^-\  for  /*,  ^L   Variable 
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Maximum  Error  of  Limit  Cycle  from  Unit  Step  Input 
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Maximum  Error  of  Limit  Cycle  from  Unit  Step  Input 
Limit  Cycle  Error  (Radians)  (^^-1  for  /*,   '   Variable 
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APPENDIX     A 

COMPUTAT/O/V      FLOW  £>/A  G-/QA/V1 

for 

CDC   /604      D/G-/TAL     COMPUTER 
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APPENDIX  A  (continued) 
EXPLANATION  OF  FLOW  DIAGFAM  BLOCK  SYMBOLS 

1.  Clear  fault  stop. 

2.  Print  last  overshoot  stored  in  LIMIT3UF  (symbol  868). 
Print  out  position  in  phase  trajectory  of  fault  stop 
(symbol  666). 

3.  End  of  file. 

^4 .  Actuate  fault  stop. 

5-   Load  zeta  from  Index  1  (b  ). 

6.  Print  last  overshoot  stored  on  LIMITBUF  (symbol  888). 
Print  position  in  phase  trajectory  where  it  was  manually 
stopped  (symbol  I96I ) . 

7.  Stop  if  Zeta  =  0. 

8.  Zf  CJh    =  a 

9.  Is   zeta  =  0.1? 

10.  Is  zeta  less  than  0.1? 

11.  Set  up  computation  time  of  0.01  seconds  for  problem  and 
set  up  for  print  at  every  0.1  seconds  vhen  zeta  is  greater 
than  0.1. 

12.  Set  up  computation  time  of  0.00^  seconds  for  problem  and  set 
up  for  print  at  every  0.1  seconds  when  zeta  is  equal  to  or 
less  than  0.1. 

13.  Set  up  limit  cycle  print-out  routine  for  20  overshoots  and 
average  the  last  8. 

1^.   Set  up  number  of  limit  cycle  print-outs  and  averaging  routine 
depending  upon  the  size  of  zeta: 
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T 

Overshoots 

Start 
Averaging 

Divide 
by  * 

Comp. 
Time 

^0.8 

6 

4 

2.0 

0=01 

0.6 

7 

5 

2,0 

0.01 

0.5 

8 

5 

3-0 

0.01 

0.4 

10 

7 

3-0 

0.01 

0.3 

14 

10 

4.0 

0.01 

0.2 

20 

12 

8.0 

0.01 

^  0.1 

17 

12 

5-0 

0.004 

*The  last  row  printed  has  this  number  printed  on  the  far  left 
if  the  problem  has  a  normal  solution.   This  number  indicates 
how  many  overshoots  were  averaged.   The  average  value  has  1.0 
subtracted  from  it  and  is  printed  in  the  far  right  of  the  same 
last  row. 

15.  Page  reject. 

16.  Clear  buffers  and  set  up  for  new  problem. 

17.  Restitution,  e.  from  Index  3  (b^). 

5 

Backlash,  &   ,    from  Index  4  (b,  ). 


4 


Jm  from  Index  5  (b^). 

5 


Jt  ~  Jm  K  Jmr\ 


X  1U1U   i.UU.CA  \J        \   U£  /  . 

ft 


IP  from  Index  6  (b^) 


18.  Have  all  values  been  cycled  through? 

19.  Set  ST0PN0W  after  this  run  to  stop  absolutely. 
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20.  Backlash,  /\,  y  =0? 

21.  Print  consecutive  number  of  this  parameter  run  and  all  the 
parameters  used: 

K      Vj^C     tyJm=D   VjJ^   *>„*     J* 

22.  Set  program  so  always  acts  as  combined  system. 

23.  Print  consecutive  number  of  this  parameter  run  and  selected 
parameters  for  linear  system! 

2k.      Solve  the  second  order  equation  by  Runge-Kutta-Gill  numerical 

integration  every  0.01  seconds  problem  time  when  zeta  is  greater 
than  0.1,  otherwise  every  0.004-  seconds.   Initial  conditions 

•  mm 

are  (SL  =  1.0,    (3)  =  Q   =  0  and   (Q  =  1.0  or  the  pre- 
viously  computed  point.   Next  point  is  computed  by  four 
iterations  of: 

ec=  <oh(eR-0c)-2/wh  ©c 

25.  Set  EXIT  from  print  routine. 

26.  Load  time,   (3c.  an(^-  St  ^n  P1*^11^  buffer. 

27.  Go  to  LIMITCYC  to  print  only  the  overshoot.   Otherwise  print 
the  phase  trajectory  points  for  each  0.1  seconds  up  to  a  maxi- 
mum of  12  cycles. 

28.  Store  0  for  ©   and  ©m 


W 
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29.  Is  problem  time  0.1  sec? 

30.  Is  problem  time  greater  than  0.1  sec? 

a 

31.  Print  the  number  of  point  computed,  time,  ©c  and   0C  . 

32.  EXIT  (automatically  set  for  desired  jump-out. ) 

33.  Store  time,  ^L  ,  ©c  and  @c  for  motor  and  load  initial 
conditions. 

3^.   Is  slope  NT  =  N  ? 

s —      f*+ y  o-'nT'  — =-  =  positive  value? 

35 •   Is  £)   negative? 

36.      Store    0C~    A  irP     0W1 

/° 
37-      Solve  @c=     -i^l&c)    for    ©c  and    ©<T 

by  Runge-Kutta-Gill. 

for   w*vi   and   ^^n*  by  Runge-Kutta-Gill 

39.  6>L  -/°  e»\  ^  o  ? 

1*0.  0L    r/O  ©  kw  "■   A     =  negative? 

4l.  Set  EXIT  from  print  routine. 

k2.  Set  EXIT  from  print  routine. 

43.  Load  time,  (^  and  Qc    , 

kk.  Load  @.   and  Q^V\    • 

k^>.  Print  number  of  point  computed  and  stored  values. 

k-6.  EXIT  (automatically  set  for  desired  jump-out). 
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kf.      Solve 
*/ 

©c 


\VA 


^o*(\+e)  +  e4/°zj^-e) 


e 


rv> 


7* 


h9. 

50. 

51. 

52. 
55. 
5t. 

55- 
56. 
57- 
58. 
59. 
6o. 
6i. 
62. 

65. 

64. 

65. 
66. 


Set  EXIT  from  print  routine. 

Is  (C)  positive? 

Set  up  to  count  one  phase  trajectory  when  @      is  negative. 

Is   ©c  >  1.0? 

Have  12  cycles  of  phase  trajectory  "been  completed? 

Arrange  for  no  counting  until  ^^  is  positive  again. 

Is  /^©^  "  ©c   negative? 

Set  up   -  (/°<3>k   *  <S>c)  . 

Is  y^©^  -  ® r         greater  than  preselected  £   ? 

Set  up  combined  system. 

Is  ©c  positive? 

Is  (O^  negative? 

Is  0L  —  |  positive? 

Is  (^_  —  )  negative? 

Is  0  c""  I  negative? 

Oc  v*^= 0? 

^c-//«<9»vy=  o? 

Is  (3C  ""/^  ^*fl  positive? 
Is  Q0   -yO  (^^positive? 
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67.  Is  Slope  NL  =  Ng  ? 

- -^ &  i°o;h  4.    Jt. 

A  J       ^  =  positive  value? 

68.  Modify  FLQATEST  to  keep  system  combined.   Takes  off  modification 
when   ©^   is  negative. 

69.  Is  Sc  —  /°  Qf^   positive? 

70.  Store  €>u  -  A   in  @ 

/° 

71.  Store  Qg   in  ©  ^ 

72.  Is  0C  -  J   negative? 

75-    Gc  -,/•©»»  =  0? 

7^.  Oc-  ^OvA=   0? 

75-   Is  @c  -  /v?  <SnoPositive? 

76.  Is  0£  -yO  ^^  positive? 

77.  Modify  FLQATEST  to  keep  system  combined.   Takes  off  modification 

o 
when   ©c  ^s  positive. 

78.  Have  1  l/2  minutes  of  real  time  elapsed  since  the  last  print-out? 

79.  Is  0c  <  1.0? 

80.  Set  up  to  print  when  in  the  fourth  quadrant. 

81.  Check  if  in  the  fourth  quadrant. 

82.  Is  the  new  &q   greater  than  the  previous  &c   • 

83.  Store  the  new  time,  &c     and-  ®C   * 

8k.      EXIT  (automatically  set  for  desired  jump-out). 

85.  Print  the  last  overshoot  stored  in  LLMITBUF  (symbol  888) 

86.  Stop  real  time  clock. 

87.  Print  position  in  phase  trajectory  where  it  was  stopped  by  the 
real  time  clock  (symbol  999 )• 
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88.  Return  to  "beginning  to  start  new  problem  if  new  parameters 
remain  to  be  evaluated.  Unconditional  stop  if  all  parameters 
have  been  evaluated. 

89.  Has  the  required  number  of  overshoots  been  evaluated  prior  to 
starting  the  averaging  routine? 

90.  Add  the  (9c  Just  evaluated  to  the  previous  sum  for  averaging. 

91.  Add  one  to  the  count  prior  to  averaging. 

92.  Clear  time  clock. 

o 
93«   Print  time,  ©c  and  ®c.  • 

9^.  Required  number  of  total  print-outs? 

95.  Clear  buffers. 

96.  Obtain  the  average  of  C£  . 

97.  Subtract  1.0. 

98.  Print  the  average  of  ($U—  1.0  and  how  many  overshoots  were 
utilized  to  obtain  the  average. 

99.  Is  @c  less  than  1.005  radians? 

100.  Clear  buffers. 

101.  Start  real  time  clock. 

102.  EXIT  (automatically  set  to  desired  jump-out). 

103.  Has  30  seconds  of  real  time  elapsed  since  the  last  print-out? 
10^4-.  Print  the  position  in  the  phase  trajectory  where  it  was  stopped 

by  the  real  time  clock  and  (3L  being  less  than  1.005  radians 
(symbol  1005) • 
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Jump  Switches: 


J  -  For  LIMITCYC ,   allows  the  maximum  overshoot  print-out  only. 

With  J  down,  a  phase  trajectory  is  printed  for  a  maximum 

of  12  cycles. 
J  -  Set  the  number  of  print-outs  depending  on  zeta.   Change 

the  number  of  overshoot  values  used  to  obtain  the  average 

value. 

J  -  Use  the  real  time  clock  for  an  automatic  recycle  to  a 
3 

new  parameter  when  have: 

a.  1  l/2  minutes  maximum  after  the  last  print-out. 

b.  30  seconds  maximum  and  @r   i-s  1*005  radians  or  less. 


Stop  Switches: 


S  -  Stops  at  the  end  of  this  parameter  run  or  stops  at  the  end 
of  the  manual  print -out. 
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APPENDIX  B 
COMPUTER  PROGRAM 


0000? 

40000 
40001 

40002 
40003 

40004 
40005 
40006 
40007 
40010 
40011 
40012 
40013 
40014 
40015 

40016 
40017 
40020 
40021 
40022 
40023 
40024 
40025 
40026 
40027 
40030 
40031 
40032 
40033 
40034 
40035 
40036 
40037 
4Q040 
'+0041 
40042 


7A  0  00070 
75  0  40701 


74 

II 

20 

50 
75 
75 
50' 

12 
20 
12 
22 
32 
32 
20 
12 
31 
22 
22 
12 
20 
20 
20 
12 
20 
-.0 
75 
50 

II 

12 
20 
12 
20 
10 

75 
50 
02 

04 
12 
20 
10 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
?0 

2£ 
20 

20 

12 
20 
20 
12 
20 
]  ? 

12 

20 
75 
50 


00100 
4OOO4 
40623 
40425 
00000 
40000 
4067  4 

00000 

40565 
40563 
40425 
40564 
40665 
40- 

40404 
40624 
40425 
633 
4Q433 
40402 
40154 

^-01T3 

40207 
40 

40400 
00000 
40437 
00000 
40652 
40616 
40662 
40617 
40672 
40620 
00010 
40362 
71000 
00000 
40427 
300 
40664 
40156 
00000 

00000 

40162 

40161 
40157 
40155 
40377 
40615 
40351 
40352 
40353 
40614 
40607 
40610 
40611 
40612 
40613- 
40621 
40326 
40335 
40622 
40  5  44 
40253 
40071 
40254 
40072 
40452 
00000 


FAULT 

START 
LAZYZETA 

MANSTART 

NEWZETA 


ZZ 


GOAROUND 


NEWINPUT 


ORG     00007 

EXF  0  00070 

SLJ  0  FAULPRIN 
ORG     40000 

EXF  0  00100 

SLJ  0  NE.V/ZFTA-1 

LDA  1  ZbiTAlNDX 

STA  0  ZETA 

ENI  1  0 
SLJ  '  0  START 

SLJ  0  MANUAL 

ENI  0  0 

LDA  0  NEWSTOP 

STA  0  STOPNOW 

LDA  0  ZETA 

AJP  0  STOPINDX 

FMU  0  TWO 

FMU  0  OMEGAN 

STA  0  A 

LDA  0  ZETAINDX+1 

FSB  0  ZETA 

AJP  0  OUIKCOMP 

AJP  2  QUIKCOMP 

LDA  0  POINT01 

STA  0  TABLES+1 

STA  0  TABLEL+1 

STA  0  TABLEM+1 

LDA  0  COUNT 10 

STA  0  TENTHSEC 

ENI  0  0 

SLJ  2  BIGZETA 

ENI  0  0 

LDA  0  STOP20 

STA  0  STOPZETA 

LDA  0  COUNT  12 

STA  0  COUNTZET 

LDA  0  EIGHT 

STA  0  DIVIDE 

ENA  0  10 

SAL  0  STOPRINT 

SLJ  4  DECOF 

ENI  0  0 

02  0  FORZEROS 

04  0  0 

LDA  0  ONE 

STA  0  UDOT 

ENA  0  0 

STA  0  0 

STA  0  THETA 

STA  0  THETADOT 

STA  0  U 

STA  0  T 

STA  0  INDEX 
STA        iXAVE 

STA  0  PRINTBUF+3 

STA  0  PRINTBUF+4 

STA  0  PRINTBUF+5 

STA  0  LIMINDEX 

STA  0  LIMITBUF+2 

STA  0  LIMITBUF+3 

STA  0  LIMITBUF+4 

STA  0  LIMITBUF+5 

STA  0  LIMITBUF+6 

LDA  0  PRINTSET 

STA  0  OK2PRINT+1 

STA  0  BOUNPRIN+4 

LDA  0  NEWLIMPT 

STA  0  L1MPRINT+1 

LDA  0  FLOTDATA  . 

STA  0  FLOATEST 

LDA  0  FLOTDATA+1 

STA  0  FLOATEST+1 

SLJ  0  INDXVALU 

ENI  0  000  00 


& 


40043 
40044 
40045 
40046 
40047 
40050 

40051 
40052 
40053 
40054 
40055 
40056 
40057 
40060 
40061 
40062 
40063 
40064 
40065 
40066 
40067 
40070 

40071 
40072 
40073 
40074 
40075 
40076 
40077 
40100 
40101 
40102 
40103 
£0104 
40105 
40106 
40107 


if  8 

40500 

LDA 

AJP 

*i  ?  0 

42£-3 

LDA 

20  0 

kl\%l\ 

00000 

STA 

7£  ft 
50  0 

SLJ 

eSi 

01  0 

40404 

01 

06  0 

00001 

06 

72  0 

40046 

RAO 

Z?8 

71000 
40412 

SW 

06  0 

00000 

06 

72  0 

40510 

RAO 

75  4 

71000 

SLJ 

01  0 

40420 

01 

06  0 

00000 

06 

10  0 

00000 

COMBINED  ENA 

20  0 

40160 

STA 

20  0 

40163 

STA 

lo    4 

60200 

SLJ 

00  0 

40153 

.  0 

00  0 

40164 

0 

75  4 

60201 

SLJ 

50  0 

OOOOO 

ENI 

75  0 

40310 

SLJ 

50  0 

00000 

ENI 

12  0 

40155 

LDA 

20  0 

40174 

STA 

20  0 

40210 

STA 

12  0 

40155 

LDA 

20  0 

40175 

STA 

33  0 

40411 

FDV 

20  0 

40211 

STA 

12  0 

40157 

LDA 

20  0 

40176 

STA 

20  0 

402  00 

STA 

33  0 

40/ 

FDV 

20  0 

40212 

STA 

20  0 

40214 

STA 

12  0 

40162 

LDA 

20  0 

40201 

STA 

33  0 

40411 

FDV 

20  0 

40215 

STA 

50  0 

00000 

ENI 

13  0 

40162 

FLOATEST  LAC 

30  0 

40664 

FAD 

32  0 

40424 

FMU 

33  0 

40161 

FDV 

31  0 

40404 

FSB 

30  0 

40423 

FAD 

22  2 

40053 

AJP 

12  0 

40161 

LDA 

22  3 

40270 

AJP 

50  0 

00000 

ENI 

10  0 

00000 

LOAD      ENA 

20  0 

40177 

STA 

20  0 

40202 

STA 

20  0 

40213 

STA 

20  0 

40216 
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60200 

SLJ 

00  0 

40172 

0 

00  0 

40203 

0 

75  4 

60201 

SLJ 

50  0 
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ENI 

75  4 

60200 

MOTOR     SLJ 

50  0 

00000 

ENI 

00  0 

40206 

0 

00  0 

40217 

0 

75  4 

60201 

SLJ 

50  0 

00000 

ENI 

13  0 

40215 

CONTEST   LAC 

32  0 

40411 

FMU 

30  0 

40201 

f;d 

22  0 

40112 

AJP 

0 
0 
4 
0 
0 
0 
0 
4 
0 
0 
0 
4 
0 
0 
0 
0 
0 
4 
0 
0 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
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0 
0 
2 
0 
3 
0 
0 
0 
0 
0 
0 
4 
0 
0 

L 
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4 
0 
0 
0 
4 
0 
0 
0 
0 
0 


DELTA 

LINEAR 

NOSINK 

FLOATEST 

DECOF 

0 

A 

U 

DECOF 

MOINERT 

0 

LINEAR+10 

DECOF 

KMOTCONS 

0 

0 

QX 

QY 

RUNGE 

TABLES 

DERIVES 

RUNGE+1 

0 

COMBPRIN 

0 

T 

TL 

TM 

UDOT 

VDOT 

RHO 

WDOT 

U 

V 

THETADL 

RHO 

W 

THETADM  ■ 

THETA 

THETAL 

RHO 

THETAM 

0 

THETA 

ONE 

OMEGANSO 

THETADOT 

A 

B 

COMBINED 

THETADOT 

MODIFY 

0 

0 

QL 

(  It 

OM 

QMM 

RUNGE 

TABLEL 

DERIVL 

RUNGE+1 

0 

RUNGE 

0 

TABLEM 

DERIVM 

RUNGE+1 

0 

THETAM 

RHO 

THETAL 

BOUNCE 
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40110 

22 

3 

40112 

AJP 

31 

0 

40410 

FSB 

40111 

22 

3 

40317 

AJP 

50 

0 

00000 

ENI 

40112 

75 

0 

40331 

BOUNCE    SLJ 

50 

0 

00006 

ENI 
LDA 

40113 

12 

0 

40200 

20 

0 

40212 

ST  A 

40114 

12 

0 

40417 

LDA 

33 

0 

40405 

FDV 

40115 

31 

0 

40407 

FSB 

32 

0 

40200 

FMU 

40116 

20 

0 

40200 

STA 

12 

0 

40664 

LDA 

40117 

30 

0 

40407 

FAD 

32 

0 

40411 

FMU 

40120 

32 

0 

40214 

FMU 

30 

0 

40200 

FAD 

40121 

33 

0 

40414 

FDV 

32 

0 

40412 

FMU 

40122 

20 

0 

40200 

STA 

20 

0 

40176 

STA 

40123   . 

13 

0 

40176 

LAC 

32 

0 

404  23 

FMU 

40124 

20 

0 

40175 

STA 

13 

0 

40214 

LAC 

40125 

32 

0 

40411 

FMU 

30 

0 

40212 

FAD 

40126 

32 

0 

40407 

FMU 

30 

0 

40200 

FAD 

40127 

33 

0 

40411 

FDV 

20 

0 

40214 

STA 

40130 

20 

0 

40212 

STA 

" 

13 

0 

40212 

LAC 

40131 

32 

0 

40422 

FMU 

20 

0 

^22^ 

STA 

40132 

13 

0 

40201 

LAC 

30 

0 

40664 

FAO 

40133 

32 

0 

40421 

FMU 

30 

0 

40211 

FAD 

40134 

20 

0 

40211 

STA 

75 

0 

40332 

SLJ 

40135 

50 

0 

00000 

COMBTEST  EMI 

12 

0 

40214 

LDA 

40136 

32 

0 

40411 

FMU 

31 

0 

40200 

Foo 

40137 

22 

3 

40276 

AJP 

50 

0 

00000 

EMI 

40140 

65 

0 

40403 

THS 

75 

0 

40301 

SLJ 

40141 

12 

0 

40200 

LDA 

20 

0 

40161 

STA 

40142 

20 

0 

40157 

STA 

13 

0 

40157 

LAC 

40143 

32 

0 

40404 

FMU 

20 

0 

40156 

STA 

40144 

12 

0 

40201 

LDA 

20 

0 

40162 

STA 

40145 

13 

0 

40162 

LAC 

30 

0 

40664 

.-.  D 

40146 

32 

0 

40424 

FMU 

30 

0 

40156 

FAD 

40147 

20 

0 

40156 

STA 

12 

0 

40174 

LDA 

*0150 

.  20 

0 

40155 

STA 

- 

12 

0 

40200 

LDA 

40151 

22 

2 

40225 

AJP 

22 

o 
5 

40255 

AJP 

40152 

75 

0 

40273 

SLJ 

50 

0 

00000 

ENI 

40153 

00 

0 

00000 

TABLES    OCT 

00 

0 

00002 

40154 

17 

7 

15075 

DEC 

34 

1 

21727 

3  BOUNCE 

0  DELTA 

3  FLOTPRIM 

0  0 

0  BOUNPRIM 

0  0 

0  THETADL 

0  W 

0  RHOSO 

0  TNERTRAT 

0  REST  I  TUT- 

0  THETADL 

0  THETADL 

0  ONE 

0  RESTITUT 

0  RHO 

0  THETADM 

0  THETADL 

0  TOTINERT 

0  MOINERT 

0  THETADL 

0  V 

0  V 

0  B 
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0  THETADM 
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0  0 
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0  THETADL 

3  COMPLMNT 
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0  TL 

0  T 

0  THETADL 

2  PLUSTEST 
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0  0 

2 

•  01 
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40171 
40172 
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20 

c 

00 
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00 
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00 

0 

88 
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0 

00 
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00 

0 
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0 
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0 

13 

0 

32 
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20 

0 

13 

0 

30 

0 

32 
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30 

0 

20 

0 

12 

0 

20 

0 

75 

0 

50 

0 

00 

0 

00 

.0 

17 

7 

34 
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_ 

00 

0 

00 

0 

00 

0 

00 

0 

00 

c 

00 

0 

00 

0 

00 

0 

00 

0 

00 

0 

00 

0 

00 

0 

CO 

0 

00 

0 

13 

0 

32 

0 

20 

0 

12 

0 

'0  0 
'5  0 

00  0 
00  0 


IV 
34 
00 
00 


00  0 
00  0 
00  0 
00  0 
00  0 
00  0 
00  0 
00  0 
00  0 


00 
00 


00  0 
13  0 
32  0 
20  0 
13  0 
32  0 
30  0 


>0000 
1000 
."-000 
00000 
00000 
00000 
00000 

ooooo 

00000 

ooooo 
ooooo 
ooooo 
ooooo 
ooooo 

40157 
40404 
40156 
40162 
40664 
40424 
40156 
40156 
40157 
49161 
60202 
OOOOO 
OOOOO 
00002 
15075 
217  27 

ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 

40176 
40423 
40175 
40176 
40200 
60202 
OOOOO 
00002 
15075 
21727 
OOOOO 
OOOOO 
OOOOO 
OOOOO 
OOOOO 

ooooo 
ooooo 

06000 

ooooo 
ooooo 
ooooo 
ooooo 
ooooo 
ooooo 

40201 
40421 
40211 
40212 
40422 
40421 


T 

UDOT 

U 

QX 

THETADOT 

THETA 

QY 

DERIVES 


TABLEL 


TM 

WDOT 

W 

QM 

THETAOM 

THETAM 

QMM 

DERIVM 


DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
LAC 

STA 
LAC 
FAD 
FMU 
FAD 
STA 
LDA 
STA 
SLJ 
EN  I 
OCT 

DEC 


TL 

DEC 

VDOT 

DEC 

V 

DEC 

OL 

DEC 

THETADL 

DEC 

THETAL 

DEC 

QLL 

DEC 

DERIVL 

TABLEM 

-LAC 
FMU 
STA 
LDA 
STA 
SLJ 
OCT 

DEC 
DEC 
DEC 
DEC 

DEC 

DEC 

DEC 

DEC 

LAC 
FMU 
STA 
LAC 
FMU 
FAD 


0 
1*0 

0 
0 
0 
0 
0 

o  u 

0  A 

0  UDOT 
0  THETA 
0  ONE 

0  OMEGANSQ 
0  UDOT 
0  UDOT 
0  U 

0  THETADOT 
0  RUNGE+2 
0  0 
2 
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0 

0 

•  0 
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0  THETADL 
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2 

|01 

0 

0 

0 

0 

0 

0 

0 

0  THETAL 

0  UDOT 
0  W 
0  D 
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40222 

30 

0 

40211 

FAD 

0 

WDOT 

\% 

0 
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402_2 

STA 

0 

WDOT 

40223 

0 

LDA 

0 

20 

0 

40214 

STA 

0 

THETADM 

40224 

75 

0 

60202 

SLJ 

0 

RUNGE+2 

50 

0 

ooooo 

ENI 

0 

0 

40225 

12 

0 

40201 

PLUSTEST  LDA 

0 

THETAL 

31 

0 

40c 

FSB 

0 

ONE 

40226 

22 

3 

40232 

AJP 

3 

PLUSTEST+5 

13 

0 

401 

LAC 

0 

THETAM 

40227 

32 

0 

40411 

FMU 

0 

RHO 

•  30 

0 

40201 

FAD 

0 

THETAL 

40230 

22 

0 

40301 

AJP 

0 

PREPLANS 

22 

2 

40241 

AJP 

2 

PLUSFORC 

40231 

75 

0 

40301 

SLJ 

0 

PREPLANS 

50 

0 

ooooo 

ENI 

0 

0 

40232 

13 

0 

40215 

LAC 

0 

THETAM 

32 

0 

40411 

FMU 

0 

RHO 

40233 

30 

0 

40201 

FAD 

0 

THETAL 

22 

0 

40235 

AJP 

0 
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40234 

22 

2 
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2 

PREPLANS 

50 

0 
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0 

0 

40235 

13 

0 

40162 

LAC 

0 

THETA 

30 

0 

40664 

40424 

FAD 

0 
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32 

0 
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0 
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23 

0 

•FDV 

0 

THETADOT 

40237 

31 

0 
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FSB 

0 

A 

30 

0 

40423 

FAD 

0 

B 

40240 

22 

2 

40053 

AJP 

2 

COMBINED 

■ 

75 

0 
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SLJ 

0 

PREPLANS 

40241 

12 

0 

40244 
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0 
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20 

0 

40071 
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0 

FLOATEST 

40242 

12 

0 

40245 

LDA 

0 
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20 

0 

40072 
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0 
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75 

0 
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0 
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50 

0 
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0 

0 
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12 

0 
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0 
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22 

2 
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2 
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75 

0 
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0 
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50 

0 
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0 

0 
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12 

0 
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0 
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22 

3 
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3 
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0 
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0 
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50 

0 
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ENI 

0 
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0 
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50 
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75 

3 

75 
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40607 
40356 
40360 
40346 

4  0 ".  0  5 

40347 
40606 
40350 
40607 
40360 
77775 
40360 
40615 
40617 
40541 
40617 
40615 
40615 
40  543 
00000 
40607 
40613 
40613 
00*300 
02000 
71000 
40605 
00001 
40544 
40614 
40614 
40616 
40555 
40616 
40614 
40604 
40607 
40573 
00000 
40607 
00000 

01000 

40360 
00000 
00000 
40605 
40606 
40607 
40613 
40620 
40664 
40610 
71000 
00000 
40605 
00000 
40022 
00000 
40000 
00000 
40022 
00000 
OOCOO 
71000 
40605" 
01570 
02000 
71000 
40346 
01747 
40563 
00000 
40574 
40552 


!:§§  8  \>mmti 

AJP  2  EXIT 

LDA  0  PRINTBUF 

STA  0  LIMITBUF 

LDA  0  PRINTBUF+1 

STA  0  Li'MITBUF  +  1 

LDA  0  PRINTBUF+2 

STA  0  LIMTTBUr+2 

SLJ  0  EXIT 

PRINTNOW  ISK  2  77775 

SLJ  0  EXIT 

LDA  0  INDEXAVE 

SUB  0  COUiMTZET 

AJP  0  /+3 

ADD  0  COUNTZET 

STA  0  INDEXAVE 

RAO  0  INDEXAVE 

SLJ  0  LIMPRINT 

EN  I  0  0 

LDA  0  LIMITBUF+2 

FAD  0  LIMITBUF+6 

STA  0  LIMITBUF+6 

ENI  0  0 

LIMPRINT  EXF  0  02000 

SLJ  4  DECOF 

01  0  LIMITBUF 

06  0  1 

RAO  0  /-l 

RAO  0  LIMINDEX 

LDA  0  LIMINDEX 

SUB  0  STOPZETA 

AJP  0  STOPAVE 

ADD  0  STOPZETA 

STA  0  LIMINDEX 

LDA  0  TOOSMALL 

FSB  0  LIMITBUF+2 

AJP  2  SMALTIME 

ENA  0  0 

STA  0  LIMITBUF+2 

STA  0  0 

EXF  0  01000 

SLJ  0  EXIT 

ENI  0  0 

STOPAVE   ENA  0  0 

STA  0  LIMITBUF 

STA  0  LIMITBUF+1 

STA  0  LIMITBUF+2 

LDA  0  LIMITBUF+6 

FDV  0  DIVIDE 

FSB  0  ONE 

STA  0  LIMITBUF+3 

SLJ  4  DECOF 

ENI  0  0 

STOPRINT   01  0  LIMITBUF 

06  0  0 

STOPNOW   SLS  1  GOAROUND 

ENI  0  0 

STOPINDX  SLS  0  START 

ENI  0  0 

NEWSTOP   SLS  . 1  GOAROUND 

ENI  0  0 

TIMESTOP  ENI  2  0 

SLJ  4  DECOF 

01  0  LIMITBUF 

04  0  1570 

EXF  0  02000  . 

SLJ  4  DECOF 

01  0  PRINTBUF 

04  0  1747 

SLJ  0  STOPNOW 

ENI  0  0 

SMALTIME  SLJ  3  /+1 

SLJ  0  LIMPRINT+7 
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03410 

00  0 

OOOOO 

TIME90SC 
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00  0 
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0*1 

63  1 

46314 

17  7 

56314 

DEC 
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